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Before a redshift z=100, about 20 million years after the big bang, the uni- 
verse was nearly smooth and homogenous. 1 After this epoch tiny fluctuations 
imprinted in the matter distribution during the initial expansion began to 
collapse via gravitational instability. The properties of these fluctuations de- 
pend on the unknown nature of dark matter, 2 4 which is one of the biggest 
challenges in present day science. 5 7 Here we present supercomputer simula- 
tions of the concordance cosmological model assuming neutralino dark matter 
and find the first objects to form are numerous earth mass dark matter halos 
about as large as the solar system. They are stable against gravitational dis- 
ruption, even within the central regions of the Milky Way, and we expect over 
10 15 to survive within the Galactic halo with one passing through the solar 
system every few thousand years. The nearest structures will be amongst the 
brightest sources for gamma-rays from particle-particle annihilation. 

The cosmological parameters of our universe and initial conditions for structure forma- 
tion have recently been measured via a combination of observations including the cosmic 
microwave background (CMB) 8 , distant supernovae 9,10 and the large scale distribution 
of galaxies. 11 Cosmologists now face the outstanding problem of understanding the origin 
of structure in the universe from its strange mix of particles and vacuum energy. 12 ' 13 

Most of the mass of the universe must be a non-baryonic particle 1 ' 14 that remains 
undetected in laboratory experiments. The leading candidate for this "dark matter" is 

the neutralino, the lightest supersymmetric particle which is predicted to solve several key 

*Present adress: Department of Astronomy Astrophysics, University of California, 1156 High Street, 
Santa Cruz CA 95064, USA 



2 Diemand et al. 

problems in the standard model for particle physics. 5 This cold dark matter (CDM) candi- 
date is not completely collisionless. It can collide with baryons thus revealing its presence 
in laboratory detectors, although the cross-section for this interaction is extremely small. 
In a cubic meter detector containing ~ 10 30 baryon particles only a few collisions per day 
are expected from the ~ 10 13 dark matter particles that flow through the experiment as 
the earth moves through the galaxy. The neutralino is its own anti-particle and it can 
self-annihilate creating a shower of new particles including gamma-rays. 5 The annihila- 
tion rate increases as the density squared therefore the central regions of the Galaxy and 
its satellites will give the strongest signal. 15-18 However the expected rate is very low - 
the flux of photons on earth is the same as we would receive from a single candle placed 
on Pluto. Numerous experiments utilising these effects are underway which may detect 
the neutralino within the next decade. 7 Furthermore, in the next few years the LHC at 
CERN will confirm or rule out the concepts of supersymmetry (SUSY). 6 

We followed the growth and subsequent gravitational collapse and virialisation of 
the first structures in the cold dark matter universe with supercomputer calculations. 
The challenge is to accurately follow the evolution of the universe on scales that are 
many orders of magnitude smaller than previously studied, whilst also capturing the 
gravitational dynamics from large scales. We use a multiscale technique 19 in order to 
achieve the desired resolution within a small average density patch of the universe which 
is nested within a hierarchy of larger and lower resolution grids of particles. 

The fluctuations are imposed on the particles using accurate calculations of the linear 
theory power spectrum for a SUSY model with a particle mass m u = 100 GeV. This 
includes collisional damping, free streaming and the transfer of fluctuations through the 
matter-radiation era of the universe. 2-4 The resulting power spectrum is close to a power 
law P(k) oc k n with n = — 3 with an exponential cut-off at 0.6 comoving parsecs which 
corresponds to a mass scale of 1O _6 M . The cutoff scale depends on the neutralino 
mass and decoupling energy. From accelerator searches we know that m v > 37 GeV and 
the cosmic matter density sets an upper limit at 500 GeV. The damping scale for the 
allowed neutralino models differ from the model we used by less than a factor of three 
in mass 2-4 and therefore structure formation is very similar in all SUSY-CDM scenarios. 
A less popular CDM candidate is the axion, it has a much smaller damping scale of 
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1O~ 13 M . For comparison we simulated the high resolution region with an axion CDM 

fluctuation spectrum on the resolved scales. Both models produce equal halo abundances 

above 5 x 1O~ 6 M , but the axion model also forms bound structures down to the smallest 

resolved scales, see Figure 3. In this letter we concentrate on the properties of the first 

structures to form in the SUSY-CDM model. 

We evolve the initial particle distribution using a parallel multi-stepping treecode, 

starting at a redshift z = 350 when the fluctuations are still linear. The high resolution 

region forms the first non-linear structures at z=60 and the entire region quickly becomes 

distorted by the complex tidal field from the surrounding overdensities. By a redshift 

z=26 the high-resolution region begins to merge into the lower resolution surroundings 

and we do not analyse the region further, however this is sufficiently late that about 5 

percent of the mass in the region has collapsed into bound dense structures (halos), see 

Figure 1. 

The first dark matter halos to collapse and virialise are smooth triaxial objects of 
mass 1O~ 6 M and half mass radii of 1CT 2 parsecs. Figure 2 shows the density profiles 
of three representative halos at z=26 which are well fitted by single power law density 
profiles p(r) oc r -7 with slopes 7 in the range from 1.5 to 2, similar to galactic halos 
shortly after their formation. 20 Note that the densities at the virial radius are about an 
order of magnitude higher than the density at 0.01r virial in a galactic halo today, which 
makes the survival of many of these halos as galactic substructure possible. The central 
resolved densities reach 10 9 times the mean background density at one percent of their 
virial radii. Unlike galaxy and cluster mass CDM halos, they do not contain substructure 
since no smaller mass halos have collapsed in the hierarchy. 

Figure 3 shows the mass function of halos. We use a friends of friends algorithm with a 
linking length set to identify the dense central regions of collapsed halos, then for each halo 
centre we recursively search for the radius r2oo that is at an overdensity of 200 times the 
cosmic mean density. The resulting halo mass function is steep dn(M) / d log M oc M _1 . 
For comparison we plot an extrapolation of the halo mass function found on much larger 
scales > 1O 7 M [ 21 ] which fits surprisingly well up to the cut-off scale of 1O~ 6 M below 
which we find no more structures. 

At these epoch the baryons are kept sufficiently warm by the CMB that they are 
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unable to cool and form visible objects such as stars or planets within such tiny systems. 13 
The dark halos may be detected via gravitational effects such as lensing or dynamical 
perturbations. Although we can not simulate the entire galactic halo at the resolution 
required to determine the survival statistics of these objects we can make some simple 
estimates of their survival and abundances. As the Galactic halo is assembled, these first 
objects merge successively into more massive systems. From scales of 1O 7 M to 1O 15 M 
the mass function of substructure is a self similar power law of slope dn{m)/dm oc m~ L9 
[ 22 ]. Extrapolating the subhalo mass function to the smallest scales gives us a total 
number of substructure halos N(M > 1CT 6 M Q ) ~ 5 x 10 15 and the expected number 
density of subhalos at the solar radius is n(R e ) ~ 500pc~ 3 , assuming that they trace 
the mass. Although this extrapolation is made to much smaller masses than simulated 
previously, the substructure within halos collapsing at z ~ 15 with masses ~ 1O 7 M fit 
the extrapolation from larger mass scales 23 even though they form from regions of the 
CDM power spectra with effective index n as —2.95. 

Can these structures survive the strong disruptive gravitational forces from the 
Galaxy? The tidal radius is simply the inner Lagrange point of the rotation of the 
two body system. For halos with power law density profiles p(r) oc r~ 2 we find 
r t = (Rvsat)/(V^V parent ) where v sat and V parent are the effective circular velocities 
(V = y/GM/r) of the satellite and main halo and R is the pericentric distance of the 
satellite. For the smallest mini-halos v sat lm/s, r 20 o = 0.01 parsecs. Therefore within 
the Galactic potential these halos could survive completely intact to about 3 kpc from 
the centre, well within the galacto-centric position of the sun. Encounters between halos 
and with stars and molecular clouds may disrupt a small fraction of these structures but 
using the impulsive heating approximation we estimate that most will survive with little 
mass loss. 

A significant fraction of the mass may lie within bound structures at our location 
within the Galaxy, lowering the available smoothly distributed matter necessary for direct 
detection experiments. The earth passes through a dark matter mini-halo every 10,000 
years, an encounter which lasts for about 50 years, therefore most of the time the earth is 
within an underdense region of dark matter. Integrating the mass function from 1O _6 M 
to 1O 1O M , normalised such that 10% of the mass is within substructure above a mass scale 
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of 10 7 M Q (as given by simulations of Galactic halos) we find that about 50 percent of the 
mass is bound to dark matter substructures. The velocity perturbation to a planetary 
orbit or satellite is very small (« lCT 6 m/s), well below the observational constraints. 
However resonant encounters and the cumulative effects of ~ 10 6 impulsive encounters 
may cause significant perturbations to some of the bodies orbiting in the Oort cloud 
surrounding the solar system. 

Compact objects in the mass range considered here could produce a microlensing signal 
in a multiply lensed quasar image, such as time varying flux differences. 24 The lensing 
object needs to be smaller than the Einstein radius 



rE = 3 - 7xl0 ' 6 fe cm • (1) 

For a 1O~ 6 M object r E — 1CT 7 pc, which is much smaller than the size of the mini-halos 
considered here, therefore it is unlikely that gravitational lensing can provide a constraint 
on their presence, either in our halo or on cosmological path lengths to distant quasars. 

Indirect detection is a more interesting possibility and several ongoing and planned 
experiments aim to detect the atmospheric Cerenkov light from gamma-rays produced by 
neutralino annihilation in the cores of dark halos halos. 7 ' 25 ~ 28 Simple scaling arguments 
can show that minihalos can have high relative luminosities in 7-rays. The absolute 7-ray 
luminosity of a dark matter halo with an NFW density profile is proportional to L oc p 2 r^, 
where r s is the scale radius of the NFW profile and p s = p(r s ). 29 The relative luminosity 
that would arrive at the detector from a halo at a distance d is then L rel oc Ld~ 2 . 

Now we compare the relative luminosity of a minihalo at a distance of 0.1 parsec (their 
expected mean separation) to the signal from the centre of the Draco dwarf galaxy: 

Wini oc ( 7x 10 6 p crit W o.005pc W 0-lpc V 2 ^ 5 (2) 

£rel,draco ^ V 1 .7 X 10 5 p C rit / \ 300 P C ) l^OOOpCyl ~ 1 ' 

where we used the typical minihalo properties from our simulations. The large abundance 
of the smallest sub clumps compensates their smaller absolute luminosity and the closest 
of them will be bright sources of 7-rays. The background flux will be enhanced by a 
boost factor of over two orders of magnitude over a smooth Galactic dark matter po- 
tential. Current indirect detection experiments such as VERITAS, 26 HESS, 27 MAGIC 28 
or CANGAROO-III 25 can probe part of the parameter space predicted by SUSY theory 
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by observing the galactic centre. However this region is dynamically complex since it 
contains numerous confusing astrophysical gamma ray sources and a supermassive black 
hole that can erase the central cusp. CDM mini-halos are potentially bright and will not 
suffer from these problems. All sky surveys could detect some nearby minihalos which 
would have a characteristic extent on the sky that is similar to that expected for a more 
distant satellite galaxy like Draco. 
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Figure 1. A zoom into one of the first objects to form in the universe. The colours show 
the density of dark matter at redshift 26. Brighter colours correspond to regions of higher 
concentrations of matter. The blue background image shows the small scale structure in the 
top cube (cube size = [3 comoving kpc] 3 ) which has a similar filamentary topology as the large 
scale structure in the CDM universe. The first red image zooms by a factor of one hundred into 
the average density high resolution region. This region was initially a cube of [60 comoving pc] 3 
resolved with 64 million particles with a gravitational softening of 10 -2 comoving parsecs and 
masses 1.2 x 1(T 10 M Q = M moon /300. The final image shows a close up of one of the individual 
dark matter halos in this region, again zooming in by a factor of one hundred so that the box 
has a physical length of 0.024 parsecs. This tiny triaxial Earth mass halo has a cuspy density 
profile and is smooth, devoid of the substructure that is found within galactic and cluster mass 
dark matter halos. Even though the index of the power spectrum is very steep on these scales, 
n ~ —3, we find that halos can collapse before merging into a larger system, rather than the 
niave expectation that all scales are collapsing simultaneously thus erasing such structures. 
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Figure 2. Radial density profiles of three typical minihalos at redshift 26. The radial distance 
is plotted in physical units and we show low concentration a/^-profiles for comparison. We use 
the mean dark matter profile infered from the highest resolution galaxy cluster simulations, 30 
i.e. (a/?7) = (1,3,1.2). The vertical dotted line indicates our force resolution and the arrows 
indicate the radii that is 200 times the background density. Across the entire range of halo 
masses from 10 -6 to 1O 1 M0, we find small concentration parameters c < 3. We do not observe 
a trend of concentration with mass, possibly because the halos all form at a similar epoch as 
expected when the power spectrum is so steep. 
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Figure 3. The abundance of collapsed and virialised dark matter halos of a given mass. The 
same region was simulated twice using different types of intial fluctuations: (A) SUSY-CDM 
with a 100 GeV neutralino (stars) and (B) an additional model with no small scale cut-off 
to the power spectrum (open circles) as might be produced by an axion dark matter candi- 
date. Densities are given in co-moving units, masses in r L M = 1.41M Q , where h = 0.71 
is the normalized present day Hubble expansion rate. Model (B) has a steep mass function 
down to the resolution limit whereas run (A) has many fewer halos below a mass of about 
5 x 10~ 6 /i M© = 3.5 x 10 _6 /i _1 Mq. (Our simulations do not probe the mass range from about 
3 x 10 /t Mq to 2 x 10 /i -1 M0.) The dashed-dotted line shows an extrapolation of the 
number density of galaxy halos (from 21 ) assuming dn(M)/d\og M oc M . The solid line is the 
function dn(M)/dlogM = 2.8 x W 9 (M / H^Mq)- 1 exp[-(M/M cutoff )- 2 / 3 ](/i- 1 Mpc)- 3 , with a 
cutoff mass M cuto ff = 5.7 x 1O- 6 /i- 1 M . The power spectrum cutoff is P{k) oc exp[— (k/kf s ] ), 
where kf s is the free streaming scale and assuming k cx M -1 / 3 motivates the exponent of —2/3 
in our fitting function. 



